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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL NOTE 3391 


FREE-FLIGHT MEASUREMENTS OP TURBULENT-BOUNDARY-LAYER 
SKIN FRICTION IN THE PRESENCE OF SEVERE AERODYNAMIC 
HEATING AT MACH NUMBERS FROM 2.8 TO ?.0 
By Simon C. Sommer and Barbara J. Short 

SUiyiMARY 


Experimental measurements of average skin friction of the turbulent 
boundary layer have been made on free -flying, hollow- cylinder models at 
Mach numbers of 2.8, 3*8^ 5*6, and 7*0, at conditions of high rates of 
heat transfer. It has been found that for these high heat-transfer con- 
ditions, the ratio of skin friction to incompressible skin friction is 
approximately 35 percent higher than zero -heat -transfer wind-tunnel data 
at Mach numbers of 2.8 and 3*8* Although no measurements of skin fric- 
tion have been made at zero -heat -transfer conditions at veiy high Mach 
numbers, the data of the present investigation indicate that this same 
trend of increasing skin-friction ratio with increasing heat-transfer 
rates will persist at Mach numbers as high as 7* 

The Rube sin and Johnson T* method of calculating skin friction for 
laminar boundary layers has been modified and compared to the data of 
this investigation and existing wind-tunnel data for conditions close to 
zero heat transfer. It has been found that values of skin-friction ratio 
computed by this method agree well with the experimental values over a 
wide range of Mach numbers and heat-transfer conditions. 


INTRODUCTION 


The present state of knowledge of* the skin friction of turbulent 
boundary layers at supersonic speeds is primarily guided by the experi- 
mental data that exist. These data are fairly complete for conditions 
close to zero heat transfer at Mach numbers up to (refs. 1 and 2). 
Unfortunately, there has been little experimental investigation of the 
effects of heat transfer and further increases in Mach number on skin 
friction. Theoretical estimates generally agree that skin friction 
increases with increasing heat transfer from the boundary layer to the 
wall, and decreases with increasing Mach number (e.g., refs. 3 through 7) ^ 
but are not generally in agreement quantitatively. Since heat-transfer 
rates will probably be large under conditions of free flight and since 
flight speeds of interest extend well beyond a Mach number of 4.5, a 


2 


NACA TN 3391 


program vas initiated in the Ames supersonic free -flight wind tunnel to 
measure skin friction of the turbulent boundary layer \mder conditions 
of large heat transfer and to extend the Mach number range for which 
skin-friction data are available. The results of this investigation are 
reported herein.^ 


SYMBOLS 


A ratio of that part of the trip drag which results in removing 

momentiJin from the houndaiy layer to the total drag, — (see 
Appendix B) , dimensionless 

total-drag coefficient, dimensionless 

CDt trip-drag coefficient, dimensionless 

coefficient of that part of the trip drag vhich results in remov- 
ing momentum from the Boundary layer (see Appendix B) , dimension- 
less 

Cjp average skin-friction coefficient, turbulent flow, dimensionless 

incompressible skin-friction coefficient, turbulent flow, dimen- 
sionless 

Cp average skin-friction coefficient, laminar flow, dimensionless 

Cm specific heat of model material, Btu/lb °F 

cp specific heat of air at constant pressure, Btu/lb °F 

H average heat-transfer coefficient, Btu/sec sq ft °F 

h wall thickness at base of model, ft 

km thermal conductivity of the'model material, Btu/sec sq ft °F/ft 
L length of run of turbulent flow, ft 

I length of model, ft 

M Mach number, dimensionless 

p static pressure, Ib/sq ft 

^Preliminary results of the present investigation have been pre- 
sented in reference 8 . 
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base pressiire, Ib/sq ft 
q dynamic pressvire, Ib/sq ft 

R Reynolds number based on model length, dimensionless 

Rl),Rc Reynolds numbers used in determining incompressible skin-friction 
coefficient, dimensionless 

Rj) Reynolds number based on pipe diameter, subsonic pipe flow, dimen- 
sionless 

r radius of model from axis to wall center, ft 

S surface area, sq ft 

Si Sutherland constant, °R 

IT 

St Stanton number, , dimensionless 

Cp^Pj^Ui^ 

T absolute temperature, °R 

Ti initial temperature of the model, °R 

t time, sec 

u velocity in the x direction of air in the boundary layer, ft/sec 

X axial distance, ft 

y radial distance, ft 

half -wall thickness, ft 


a 


n 


6 

e 

U 

P 


thermal diffusivity of the model material, — sq ft/sec 

^mPm 

yiH 

positive roots of p tan 3 = -r — (values tabulated in Appendix IV , 
ref. 2h) , dimensionless ™ 

boundary-layer thickness, ft 

boundary-layer momentum thickness, / (l ~ ft 

Pi ^1 \ Ui y 

coefficient of viscosity, lb sec/sq ft 
density of air, Ib/cu ft 
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Pjjj density of the model material, Ib/cu ft 

Subscripts 

Except where otherwise defined, the following subscripts apply: 

0 free -stream conditions 

1 conditions at the outer edge of boundary layer 

w conditions at wall 

Superscript 

’ conditions at which incompressible flow relations must be evalu- 

ated in order to represent compressible flow 

EQUIIMENT AND TEST CONDITIONS 


Skin friction was obtained from measurements of the total drag of 
spin-stabilized thin-walled tubes of the type shown in figure 1. Test 
and tare models, identical except for length, were gun-launched under 
the same conditions, and total-drag coefficients were computed from 
deceleration data. Deceleration of a model was computed from its time- 
distance history which was recorded by a chronograph and shadowgraphs 
(ref. 9)* The difference between the total drag of a test model and the 
total drag of a tare model is, except for small connections, a measure 
of the average skin-friction drag of the added length of the test model. 
This tare -drag method of obtaining skin friction and this hollow-cylinder 
model configuration were chosen because only small corrections were 
required for the evaluation of skin friction. In addition, direct cor- 
relation could be made with flat-plate results inasmuch as the flow 
closely resembled two-dimensional flow (bounda 2 ny-layer thicknesses were 
small compared to the radius of the cylinder) . 


Models and Model Launching 


The models were made of 75 S-T aluminiim, with 1.44 inches outer 
diameter and 0.030-inch-thick walls. The outer and inner surfaces were 
polished with successively finer polishing papers, the last being 4/0 
polishing paper. The finish of some typical models observed with an 
interferometer (ref. 10 ) showed the magnitude of the peak to valley 
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roughness to he approximately 20 microinches; however, the root-mean- 
square height of the surface irregularities would he considerably smaller. 
Three nose contOTirs, a double wedge with a half -angle of 10°, a double 
wedge with a half -angle of 15°, and a circular-arc profile with a tangent 
lialf -angle at the tip of 20°, each having a leading -edge thickness less 
than 0.001 inch were used. Boundary -layer trips (fig. 2), used through- 
out to promote turbulent flow over the surface of the models, will be 
discussed later. 

The test models were 2.0 and 2.5 inches long. Tare models were 0.5 
inch long. Longer test models were desired to produce a higher percentage 
of skin friction to total drag^ but the model lengths were limited by two 
factors. At the lower Mach numbers, the model length was limited to 2.0 
inches to prevent the shock wave from the leading edge from impinging on 
the inside of the cylinder wall and causing interference with the boundary- 
layer flow. This limitation was of no consequence at the higher Mach num- 
ber because of the smaller shock-wave inclination angle, but models over 
2.5 inches long failed structurally due to extremely high acceleration 
loads encountered in the gun. 

The models were launched from a standard 37-™i rifled gun, the twist 
of the rifling being one tiorn in three feet corresponding to a helix 
angle of approximately 7-5°. The models, which were approximately 0.02 
inch smaller than the land diameter of the gun, were protected from the 
rifling by plastic film wrappers which broke away from the models as 
they emerged from the gun. A rifled aluminuim disc, used to produce spin, 
was followed by a Neoprene seal which prevented powder gas leakage . A 
test-model assembly is illustrated in figure 3. Muzzle velocities of 
3200 and 4^f00 feet per second were obtained by varying the powder charges 
and resulted in peak accelerations of 130,000 and 250,000 g's on the 
models. The structural failure mentioned previously was observed on 
models 2.5 inches long at accelerations above 300,000 g*s. 


Test Conditions 


Tests were conducted at nominal Mach numbers of 2.8 and 3.9 by fir- 
ing through still air at one atmosphere pressure where the free -stream 
static temperatiire was equal to ambient temperature. A nominal Mach num- 
ber of 7.2 was obtained by firing upstream through a Mach number 2 air 
stream where the free -stream^ static temperature was about 56 percent of 
the ambient temperature. The free-stream Reynolds number range for the 
tests was from 3^0® to 9^0®. The actual free-stream Mach numbers and 
Reynolds numbers of these tests are listed in columns (l) and (2) of 
table I. It is shown in columns (3) and (4) that the Mach numbers and 
Reynolds numbers at the outer edge of the boundary layer, denoted by the 
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Free-stream conditions, as used herein, are the properties of air 
corresponding to a stationary model in an airstream at the nominal Mach 
number of the test. 
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subscript 1, are lower than the free-stream values. These changes were 
due to the flow over the nose profiles and were calculated hy two- 
dimensional shock -expans ion theory which entails simply the calciilation 
of flow at the nose tip with the oblique -shock-wave equations and flow 
downstream of the nose tip with the Prandtl -Meyer expansion equations. 

It has been shown in reference 11 that this method is applicable over 
almost the entire region of completely supersonic flow. The calculated 
static pressure at the beginning of the cylinder was very nearly equal 
to free-stream static pressure. It was assumed that the static pressure 
had returned to free-stream value at the surface trailing edge. In order 
to represent conditions over the entire cylinder, the mean static pres- 
sure was used to calculate the values of Mi and Ri listed in table I. 

The high-heat -transfer conditions of the tests are implied by the 
difference between recovery temperatures which were of the order of 
1300° to 3200° R and the initial wall temperature of approximately 530° R 
inasmuch as the wall -temperature rise during the short flight time of 0.01 
second was only about 15° to 45°. The heat-transfer conditions are indi- 
cated by the difference between columns (9) > Tr/Ti, and (8 ), Tw/Ti, of 
table I. The values listed in column (8) were estimated theoretically by 
the method explained in Appendix A. The values listed in col\min (9) were 
calculated using a recovery factor of O.89. 


Boundary-Layer Trips 


Boundary-layer trips were applied to the inside and outside surfaces 
of each model to insure a turbulent boundary layer over the model sur- 
faces. The types of trips used are illustrated in figure 2 and consisted 
of threads and raised wedges continuous around the circumference of the 
models. These two-dimensional boundary-layer trips, although not as 
effective as three-dimensional trips, could be more easily reproduced 
and machined more accurately than any type of three-dimensional roughness. 
The trip strength for each test condition was varied until the least dis- 
turbance which consistently caused turbulence to occur on or near the 
trailing edge of the roughened region was found; that is, until no laminar 
flow was observed in the shadowgraphs behind the trailing edge of the 
trip. For illustration, figure k shows a comparison between the type of 
flow observed when a trip of 0.003-inch-deep threads was used (fig. 4(a)) 
and when a trip of 0.001-inch-deep threads was used (fig. 4 (b)) at a Mach 
number of 3.9. Careful examination of figure 4(a) reveals Mach waves 
produced by the turbulent boundary layer as far forward as the shock wave 
from the trip, indicated by the leader. The Mach waves do not appear in 
figure 4(b), and moreover, with closer observation it can be seen that 
the boundary layer over the surface of the model is laminar and turbu- 
lence probably started in the annular wake; therefore, the boundary layer 
between the leading edge and the 0.003-inch-deep trip was necessarily 
laminar and transition from laminar to turbulent flow occurred somewhere 
on the trip. It is interesting to note the shock patterns associated 
with the inside and outside flows and that there is no shock impingement 
on the boundary layer. 
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A model with 0 . 003 -inch threads, the type 'vdiich produced turbulent 
flow at Mo = 3-9 (fig. i<-(a)), was tested at Mq = 7-2 and Rq = 7x10®, 
and it was found that the boundary layer was completely laminar. This 
was particularly surprising since the former test was conducted in still 
and therefore disturbance-free air, whereas the latter test was run in 
the presence of air-stream turbulence and shock waves. When the strength 
of the trip was increased from 0.00 3 -inch to 0.010 -inch-deep threads, 
there was partial laminar flow over the surface of the model. Even a 
single annular ring of the type shown in figure 2(c), raised 0.010 inch 
above the surface w’as not satisfactory, but it was found that the double 
annular ring O.OO8 inch high (fig. 2 (c)) consistently produced tTirbiolence 
in the vicinity of the trip. It was also found that 0.006-inch-deep 
threads on a circular-arc nose profile (fig. 2(b)) would consistently 
produce turbulence in the vicinity of the trip. The effect on Mach num- 
ber of using the circular-arc nose profile was to reduce Mi to 5.6. 


DATA REDUCTION 


For a more general comparison of these data with results of theory 
and experiment, it was necessary that the ratio of skin friction to 
incompressible skin friction, C-p/C-p^, be presented. It was therefore 
required that both the skin friction and the corresponding incompressible 
skin friction be determined. 


Determination of Skin-Friction Coefficient 


Skin-friction drag was obtained by subtracting the total drag of a 
tare model from the total drag of a test model. The measured total drags 
had to be adjusted because of small variations in model geometry and test 
conditions between test and tare models . These adjustments to total 
drags were made by determining the differences in individual drag com- 
ponents of a test-tare combination. These differences were obtained from 
available theoretical and experimental information and are detailed in 
the following paragraphs. The total effect of making the aforementioned 
adjustments to the total-drag data changed the skin-friction results, 
that is, the test minus tare drag, by only 5 percent and noticeably 
decreased the spread of the total-drag data. 

Base drag .- Because of the difference in the lengths of the test 
and tare models, the boundary-layer thickness at the trailing edge was 
different, and accordingly the base drag was different; consequently, 
the base drag had to be adjusted to account for this difference in 
boundary-layer thickness. The data of Chapman, Wimbrow, and Kester 
(ref. 12) for bliint trailing-edge wings are reproduced in figure ^{a.) , 
These data, extrapolated to a Mach number of 8, are cross plotte(i in 
figure 5 (l>) . Although the extrapolation appears somewhat arbitrary, it 
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is to "be emphasized that at the very high Mach numhers, the hase drag is 
only a very small percentage of the total drag, so that fairly large 
errors could he accepted without introducing significant errors in skin 
friction. If the extrapolation is in error hy as much as 20 percent at 
the high Mach mmihers, this would result in a change of the final skin 
friction of the order of 1 percent. The values of hase drag are better 
known at the lower Mach numbers as little or no extrapolation is neces- 
sary in figure 5- An error in the hase -drag corrections of 10 percent 
at these Mach numbers would result in a change of the final skin friction 
of the order of only 1.5 percent. 

Wave drag . - Wave drag consisted primarily of two components, drag 
due to leading-edge thickness and drag due to nose profile. 

Leading-edge thickness had to he measured very carefully since the 
drag associated with leading-edge thickness ranged from 8 percent to 
l6 percent of the skin-friction drag. A metallurgical microscope at 
X5OO was used to measure the leading edges since the thicknesses ranged 
from 0.0003 to 0.0010 inch from model to model. The leading edges were 
ground flat to insure uniformity and to improve the accuracy with which 
they could he measiired. The drag associated with leading-edge thickness 
was calculated with the aid of tables in reference I3 and the assumption 
that the pressure on the frontal area was the arithmetic mean of the 
total pressure and 'static pressure behind a normal shock wave.^ The maxi- 
mum error in the measurements of leading-edge -thickness could cause 
25-percent change in the corrections due to differences in leading-edge 
drag which would result in 1. 5 -percent change in the final skin friction. 

Small differences occurred in wedge angle and circular-arc profile. 

The wave drag due to the nose profile of each model was calculated hy the 
shock-expansion method, and corrections were applied for the geometric 
differences that were measured. The maximum error in the measurements 
of geometry could cause 10-percent change in these corrections which 
would result in negligible change (0.4 percent) in the final skin friction. 


^This method of calculating the leading-edge drag was confirmed 
experimentally. A 0. 5 -inch -long hollow-cylinder model having a double 
wedge with a 10^ half -angle, 0.125-inch-thick wall, and leading-edge 
thickness of 0.064 inch was tested at a Mach number of 3 ^'^d Reynolds 
number of 10®. The svna of the estimated base drag, laminar skin friction, 
and pressure drag on the 10^ wedge were only 20 percent of the measured 
total drag. The remaining 80 percent of the total drag was the drag due 
to the leading-edge bluntness. This component of drag corresponded to an 
average pressure on the frontal area which was the mean of the total 
pressure and static pressure behind a normal shock wave. If either the 
total pressure or the static pressure behind a normal shock wave were 
used to evaluate leading-edge drag, the estimated total drag would have 
differed from the measured value by ±10 percent; whereas, if the mean 
pressure were used, the estimated total drag would have been only 1 per- 
cent different than the measured value. 
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Drag due to angle of attack .- The attitude history of each model in 
flight was detennined from the shadowgraphs, and the angle of attack was 
plotted along the flight path. The mean-square angle of attack deter- 
mined from these plots was multiplied hy the lift-curve slope to give 
the increment of drag necessary to adjust the total drag to that of a 
model at zero angle of attack. The maximum increment of drag due to 
lift was so small that a 20 -percent change in the correction applied 
would result in negligible change (O.5 percent) in the final skin fric- 
tion. It is interesting to note that no visible differences in boundary- 
layer thickness around the model were ever observed in the shadowgraphs 
at any angle of attack which occurred in the test program. 

Drag due to boundary-layer trips . - Differences in the height of the 
boundaiy-layer trips were of the order of O.OOO3 inch. Corrections due 
to these differences could not be ascertained quantitatively. Although 
it was assumed that the drag due to the trip on a tare model was equal 
to the drag of the trip on a test model and no corrections were applied, 
the differences in drag due to the geometric differences of the trips 
may have been significant. The effect on skin friction of disregarding 
the possible differences in trip drag will be discussed in the section 
on test results. 


Detennination of Incompressible Skin-Friction Coefficient 


The effective point of origin of the turbulent boundary layer must 
be known in order to determine the Reynolds niombers associated with 
experimental skin friction and to calculate the corresponding incompres- 
sible skin-friction coefficient. The position of the effective turbulent 
origin can be calculated by assuming the momentum thickness of the tur- 
bulent boundary layer equal to that of the laminar boundary layer at the 
transition point. When the boundaiy-layer trips are small, this method 
is believed to be valid and was applied as described in Appendix B. It 
was found, however, that very large trips were necessary at the higher 
Mach numbers to insure the occurrence of turbulence in the vicinity of 
the trip. In this case the assumption of equal laminar and turbulent 
momentum thicknesses at the transition point is unjustified. The effects 
of such large trips were investigated, and the methods used to treat 
these cases are reported in the Results and Discussion section. 

The incompressible skin-friction coefficient was determined from 
the Karman-Schoenherr equation by the following method. The effective 
turbulent origin was calculated with the assumption that the momentum 
thicknesses of the laminar and turbulent boundary layers were equal at 
the transition point which was estimated to be at the midlength of the 
trip (see Appendix B) . Knowledge of the position of the effective tur- 
bulent origin was applied to determine the Reynolds numbers required in 
computing the incompressible skin-friction coefficient. The experimentally 
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determined skin-friction coefficient was based on tbe area included 
between tbe indicated points b and c of the sketch; thus, the incom- 
pressible skin-friction coefficient based on the same area was required. 


/ 


Effective 


turbulent origin 

•Base of tore model 


Bose of test model. 



This was determined by subtracting the incompressible skin-friction drag 
of the area indicated between o and b from the drag of the area indi- 
cated between o and c . 

(cpi') (S)q(, - (^^ob " ^^^bc 

\ ^oc ''ob ^ ^bc 

Now ( C-ff ) and ( Ct? ) can be calculated from the Karman-Schoenherr 

V-i/ob V i4c 

equation if the Reynolds number based on the length ob, Rb and the 
Reynolds number based on the length oc , Rc are known . The notation 
Rb, Rc "Will te used hereafter to designate the Reynolds numbers used in 
determining the incompressible skin-friction coefficient. 


Accuracy 


A major effort was expended to control model geometry to insure 
small scatter in the results. In spite of the best efforts to control 
model geometry, large scatter was still present. The spread of total- 
drag results for test and tare models was approximately 5 percent, but 
the inaccuracy of each individual drag result due to errors in time- 
distance measurements and air-stream calibration was only about 2 percent, 
so it is believed that differences in drag did indeed occur due to some 
unknown cause. This scatter was greatly magnified when test and tare 
models were combined to obtain skin-friction drag, for the percentage of 
skin-friction drag to total drag varied from 30 to 60 percent, depending 
on the test condition. Because of this persistent scatter it was manda- 
tory that a large number of results be obtained at each test condition 
so that the effect of random scatter due to uncontrolled features of 
model geometry could be minimized by averaging results. It is believed 
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that the scatter that remained after all corrections were made was due 
to the presence of the houndary-layer trips, because these trips could 
not he reproduced, in all cases, from model to model with the desired 
precision. 

The table below shows the scatter in skin-friction results at each 
test condition, the number of test results (i.e., with approximately an 
equal number of test and tare models, the indiscriminate combination of 
all test and tare models at each test condition) , and the symbol that 
identifies the model geometry in table I . 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

Ml 

(Ri/x)l0‘® 

Number of 
results 

RMS 

deviation 
from mean, 
percent 

Maximum 
deviation from 
mean, percent 

Symbol 

(see table l) 

2.81 

1.58 

12 

2.4 

+h.0, to -4.2 

0 

3.82 

2.13 

12 

1.1 

+1.8, to -1.8 

□ 

5.63 

1.83 

2 h 

10.0 

+21.2, to -13.3 

0 

6.90 

2.08 

20 

8.2 

+12.1, to -16.5 

□ 

7.00 

3.1^ 

6 

9.0 

+l4.9^ to -11.2 

■ 

3.78 

2.13 

2 k 

7.1 

+13.4, to -I3.3 

> 

3.67 

1.99 

3 

2.9 

+3.0, to -3.5 

< 


Colvimn (i+) in the above table shows the root-mean-square deviation from 
the arithmetic mean of the data and, statistically, best represents the 
uncertainty in the skin-friction data. Column (5) shows the maximum 
deviation from the arithmetic mean and is not nearly representative of 
the reliability of the measurements. These were included in the table 
to show the full range of the data even though, in some cases, the result 
for one combination of test and tare model was approximately 5 percent 
beyond the range of all other results at that condition. 

The three conditions for which the scatter was the least were those 
where 0.003-inch threads were used as boundary-layer trips. The condi- 
tions for which the scatter was appreciable were those where large 
boundary -layer trips were used. The skin-friction results could be 
appreciably altered by small changes in trip geometiqr that could not be 
detected (on the inner surface of the model, for example). Row (3)^ 
where the largest deviations are shown, presented the greatest problem. 
The 0.006-inch threads used for boundary-layer trips were cut on the 
circular-arc profile, so that small errors in machining would show up as 
changes in pressTire drag of the trip and of the entire contour as well. 
These changes in geometry were not accounted for and may possibly explain 
the magnitude of the deviations. 
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Because of the large scatter in the initial results at Mi = 6.9^ 
an independent test with the same model geometry at a slightly higher 
Reynolds number was run at Mi = 7-0. The agreement between the mean 
results at the two different Reynolds numbers was better than expected, 
and the difference was well within the root -mean -square deviation of the 
data. 


The more than 3 OO shadowgraphs of the data models showed turbulent 
flow over the external surfaces almost as far forward as the boundary- 
layer trips with no intermittently turbulent and laminar boundary layers. 
Unfortunately, the boundary layers on the inner surfaces of the cylinders 
could not be examined, but it is unlikely that these boundary layers 
would be different than on the outer surfaces because the machining pro- 
cess was the same for both surfaces. The flow over the boundary-layer 
trips and about 0.1 inch downstream of the trips could not be observed 
in the shadowgraphs and possibly was intermittently laminar and turbu- 
lent. It must be emphasized, however, that this region which could not 
be observed was restricted to a small section of the tare portion of the 
test models; therefore, the errors which could be introduced, although 
not considered, would be small. 


TEST RESULTS AMD DISCUSSION 
Skin-Friction Ratio 


The experimentally determined values of skin-friction coefficient 
in the form C^/C-p^, hereafter called skin-friction ratio, are plotted 

as a function of Mach number. Mi, in figure 6. Values of skin-friction 
coefficient, Cp, are mean values of many measurements at each condition, 
and values of incompressible skin-friction coefficient, Cp^, were obtained 

from the Karman-Schoenherr equation (ref. lU) . These data are uncorrected 
for thickening of the turbulent boundary layer due to the boundary -layer 
trip. The turbulent origin was determined by assuming that the momentum 
thickness of the turbulent boundary layer was equal to that of the laminar 
boundary layer at the transition point (midlength of trip) . The values 
of Cp are unaffected by this assumption, since they were determined 
experimentally. With the assimption of no initial thickness due to the 
boundary-layer trip, or so-called natural transition, the resulting values 
of skin-friction ratio, Cp/Cpj^, can be considered as lower limits. Any 
initial thickness of the boundary layer would mean a higher effective 
Reynolds number, therefore a lower value of Cpj^ and, consequently, a 
higher value of skin-friction ratio, Cp/Cp^. 

For comparison, the force data of Chapman and Kester (ref. l) and 
Coles (ref. 2) at zero-heat -transfer conditions are also included in the 
figure. The curve drawn through these points is the mean zero -heat - 
transfer curve suggested in reference 1. A comparison of the zero -heat- 
transfer data and the uncorrected data of the present investigation 
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indicate a strong dependence of skin-friction ratio on wall -temperature 
ratio, Tv/Ti, at a given Macti mmiber. 

Effect of large boundary-layer trips on skin-friction ratio . - Since 
it was found that large "boundary-layer trips were necessary at tfech. num- 
bers of 5*6 and 7*0 to promote ttirlulent flow over the surfaces of the 
models, it was required that the effect of these large trips on "both skin 
friction and the corresponding incompressible skin friction be investi- 
gated. 

To investigate experimentally the effect of using a large trip, the 
skin friction was measured on a model with a boundary-layer trip of 
0.006-inch-deep threads and compared to results obtained from a model 
with 0.003-inch -deep threads at a Mach number of 3*8 where the experi- 
mental scatter of skin friction was small. It was found that skin fric- 
tion was 10 percent lower on the 0.006-inch threaded model than on the 
0.003-inch threaded model. The results are plotted in figure 6 as skin- 
friction ratio. The low value of skin friction was not particularly 
surprising since the apparent boundary-layer thickness* on the 0.006-inch 
threaded model, as deteimined from the shadowgraphs, was meastirably 
thicker than that on the 0.003-inch threaded model at corresponding 
stations along the cylinder. Since the boundary-layer thicknesses were 
different, the Reynolds numbers were different; consequently, the values 
of incompressible skin friction were different. In order to determine 
the effect of the large trip on incompressible skin friction, the follow- 
ing procedure was used. It appeared that there would be a relationship 
between the trip drag and the thickening of the boundary layer. The trip 
drag must result in momentum being removed from the air flowing around 
the model , and part of this momentum will be removed from the boundary- 
layer air. The amount of increase of the momentum thickness of the tur- 
bulent boundary layer is dependent upon the amount of this momentum 
change confined to the boundary layer. It was necessary, therefore, to 
make some assumption regarding the amount of momentum removed from the 
boundary layer as a consequence of the trip drag. As a limiting case, 
it was assumed that all of the momentum change due to the drag of the 
trip was confined to the boundary layer. The length of run of turbulent 
flow necessary to produce the increased momentum thickness was then cal- 
culated as explained in Appendix B. The effective turbulent origin was 
established by this length of run, and a new value of incompressible 
skin friction was determined. Skin-friction ratio for the 0.006-inch 
threaded model, corrected in this way for thickening effect of the 
boundary-layer trip, is plotted in figure 7. The data uncorrected for 
boundary-layer thickening are also shown in the figure. It can be seen 
that when only the incompressible skin friction is modified, the results 
of the 0.003-inch and 0.006-inch threaded models are in reasonable agree- 
ment (within the scatter of the experimental results of the 0.006-inch 
threaded model) . The data at Mach numbers of 5*6 and 7*0 were corrected 


*The actual boundary-layer thicknesses could not be determined due 
to diffraction and refraction effects . 
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for boundary-layer thickening effect and are also shown in figure 7 - 
The correction at Mach number of 5-6 was small ^ but the data at Mach 
ntimber 7 were raised by about 10 percent, due to the correction in 
incompressible skin friction. 

Since this momentum method involves the assimiption as to the amount 
of momentum change due to the drag of the trip being confined to the 
boundary layer, it was desired to obtain further evidence as to the 
reliability of the assumption of total moment-urn change being confined to 
the boundary layer. An independent method of determining the effect of 
large boundary-layer trips was used at a Mach number of 3.8 -vdiere appar- 
ent bo-undary-layer thicknesses could be measured from the shadowgraphs. 
The apparent thicknesses of the turbulent boundary layers of the 0 . 003 - 
inch and 0.006-inch threaded models were measured and plotted as a func- 
tion of position from the leading edge. It -was then ass-umed that -when 
the apparent thicknesses of the bo-undary layers were the same, the actual 
thicknesses were the same . The additional length of run of t-urbulent 
flow necessary to produce the boundary layer observed on the model -with 
0.006-inch threads -was found by moving the x axis until the apparent 
boundary-layer thicknesses coincided. The effective turbtilent origin 
■was then established by this amount of movement of the x axis, and a 
new value of incompressible skin friction was determined. This value of 
incompressible skin friction -was in excellent agreement with the value 
determined by ass-uming total moment-um change due to the drag of the trip 
being confined to the boundary layer. Unfortunately, no check could be 
made on the corrections used at a Mach nvimber of 7 so the final skin- 
friction res-ults are presented in figxire 7 s-s bars and are tabiilated in 
table II as a range of possible values, depending upon the amo\int of 
correction to incompressible skin friction. The agreement in skin- 
friction ratio between the 0.003-inch and 0.006-inch threaded models, 
where the turbulent origins were determined assviming normal **fully^ 
developed" t\irb-ulent boundary layers, indicates that the characteristics 
of the turbulent boundary layer over the portion of the model for which 
skin friction -was measured were not significantly affected by the pres- 
ence of the boundary-layer trips. 

Also sho-wn in figure 7 nre the zero-heat-transfer data replotted 
from figure 6. At Mach numbers of 2.8 and 3 - 8 ^ "the skin-friction ratios 
at Tv/Ti = 1.03 and I.05 are approximately 35 percent higher than the 
zero -heat -transfer data. Although no measurements of skin friction have 
been made at zero -heat -transfer conditions at very high Mach numbers, 
the data of figure 7 indicate that the same trend of increasing skin- 
friction ratio -with decreasing wall-temperature ratio will persist to 
a Mach nmber of 7 - 

Effect of spin on the present results . - Because figure 7 shows an 
appreciable effect of wall -temperature ratio on skin-friction ratio , ^ and 
since the models of the present investigation were spin stabilized, it 
■was essential to determine to what extent these res\ilts had been influ- 
enced by model spin. The possibility that rotation of the models may 
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have ^ influenced the boundary layer, thus changing skin -friction drag, 
Resent in spite of the fact that the circumferential velocity ras 
small compared to the forward velocity. For this reason, an investiga- 
tion was made by firing aerodynamically stable models of the type shown 
in figure 2(e) . It was felt that it was not necessary to repeat the 
entire investigation, and that one test condition would suffice to demon- 
strate the effect of spin on the present results. The test condition 
chosen was Mq = 3-9^ where experimental scatter was least and trip drag 
negligible . — 

Since it was necessary to move the center of gravity of the model 
as far forward ,as possible to stabilize the model in flight, it was 
necessary to use a leading-edge configuration that was fairly blunt; 
hence, the 15° half -angle was used in lieu of the 10° half -angle of 
figure 2(a). In addition, since a tare model could not be stabilized 
in flight, the skin-friction drag was obtained by estimating all of the 
other drag components: laminar skin-friction drag, base drag, drag due 

to leading -edge bluntness, wave drag, drag due to angle of attack, and 
drag due to boxondary-layer trip. It was found that at the two other 
conditions (Mx = 2.8l and = 3*82) where 0.003-inch thread trips were 
used and the trip drag was negligible, the skin-friction drag could be 
estimated to within 2 percent of the experimental values if a1 l other 
drag components were estimated and subtracted from the measured total 
drag, thereby demonstrating the reliability of this procedure. 

The result of this investigation is tabulated in row (7) of table II 
and is shown in figure 7. It can be concluded that for the spin rates 
used, the effect of spin on the present resiilts is small, probably within 
the scatter of the experiment. 


The effect of wall-temperat\ire ratio on skin-friction ratio .- The 
results of figure 7 confirm^ qualitatively the conclusions of many analy- 
ses; namely, those of von Karman, Monaghan, Tucker, Van Driest, and 
Clemmow, among many (refs. 3 through J) , that skin-friction ratio 
increases with decreasing wall -temperature ratio, at a given Mach number 
and Reynolds number. These analyses differ widely only in the magnitude 
of this increase, as is shown in figure 8, where skin-friction ratio, 
Cp/Cp^, is plotted as a function of wall-temperature ratio, T^/Ti. It 
is interesting to note that three of these theories (those of von Karman, 
Monaghan, and Tucker), predict no effect of Mach number on skin-friction 
ratio at a given wall -temperature ratio, and that wall -temperature ratio 
is the controlling parameter. Also indicated in the figure is the point 
at which wall -temperature ratio is equal to recovery-temperature ratio, 
^w/Ti = Ti*/Ti, since the difference between the recovery-temperature 
ratio and the wall -temperature ratio is a meas\ire of the rate of heat 
transfer through the boundary layer. For consistency and simplicity, 
values of skin-friction ratio from these theories were calculated with 
the assumption that (|i) ~ (t)®*"^®, at a Reynolds number R^ = lOj. Incom- 
pressible skin-friction coefficients were calculated from the Karman- 
Schoenherr equation. 
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The experimental data of figure 7 were included in figure 8 to check 
the reliability of any one theory to predict skin-friction ratio over a 
wide range of Mach numbers and heat-transfer conditions. The theories 
of Van Driest and Clemmow agree well with experimental data at wall- 
temperature ratios near unity. As wall temperature approaches recovery 
temperature, these theories overestimate the skin-friction ratio by as 
much as 20 percent. Although the theories of Tucker and von Karman were 
devised for application at zero -heat -transfer conditions (Tv/Ti = Tr/Ti) , 
it is interesting to evaluate Cp/Cpj^ at conditions of large heat trans- 
fer and compare these results with experiment. At a wall -temperature 
ratio of unity, the theories of Tucker, von Karman and Monaghan predict 
a skin-friction ratio of unity, which greatly overestimates skin-friction 
ratio. At recovery temperatures where experimental data are available. 
Tucker and Monaghan predict skin-friction ratio reasonably well. 


The T' Method For Evaluating Skin Friction 


Since no one theory for turbulent flow adequately predicts the 
effects of both wall -temperature ratio and Mach number on skin-friction 
ratio, the present authors used the T’ method of Rubesin and Johnson 
(ref. 15) which was developed for laminar flow. Fischer and Norris 
(ref. 16) in 1949 applied the results of Rubesin and Johnson to corre- 
late heat-transfer data for turbulent flow. Although there was no appar- 
ent reason for choosing this laminar-flow method for turbulent flow, 
Fischer and Norris found that their data when evaluated on both the T* 
and Tv bases correlated better than when evaluated on a Ti basis. 

They indicated that it was inconclusive as to whether the T’ method of 
reference 15 should be used for turbulent flow. Young and Janssen (ref. 
17) in 1952 applied the T’ method of Rubesin and Johnson to evaluate 
skin-friction ratio for turbulent flow. They compared skin-friction 
ratio determined by this method with some zero -heat -transfer data over 
the limited Mach nimiber range from 1.5 to 2.5. The agreement was good. 

With the more complete data now available at zero-heat-transfer con- 
ditions and with the results of the present investigation at large rates 
of heat transfer, it is the purpose of this section to demonstrate the 
reliability of the T* method in predicting the effect of heat transfer 
as well as Mach number on skin-friction ratio for turbulent flow. 

The T’ method of Rubesin and Johnson consisted of finding a tem- 
perature, T’, at which the density and viscosity for compressible flow 
must be evaluated if incompressible flow relations for zero heat transfer 
are to apply. The theoretical results of Crocco-Conforto for laminar 
flow on a flat plate were used to find this reference tempera tirre , T* . 
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The following expression was derived: 


— = 1 + 0.032 Mi^ + 0.58 


( 



) 


( 1 ) 


The present authors applied equation (l) and the Karman-Schoenherr incom- 
pressihle skin-friction equation for turbulent flow to determine whether 
equation (l) wo\ild predict skin-friction ratio as a function of both Mach 
number and wall-temperature ratio. Curves of skin-friction ratio as a 
fimction of Mach ntunber, calculated from equation (l), are plotted in 
figure 9 , and are compared to zero-heat-transfer data and the primary 
data at Tv/Ti = I .03 and Tv/Ti = I .05 from figure ?• Although equa- 
tion (1) underestimates the skin-friction ratio at zero heat transfer 
above a Mach number of 2.5, the shape of the zero -heat -transfer curve is 
remarkably similar to the trend of the experimental data. Equation (l) 
also predicts, reasonably well, the effect of large rates of heat trans- 
fer on skin-friction ratio. 

Encoiiraged by these resiilts, the authors found new coefficients for 
equation (l) using the data at a Mach number of 3»82 from the present 
experiment and from the zero -heat -transfer curve. The following expres- 
sion was obtained: 


It is interesting to note that the constants of equation (2) evaluated 
for turbiilent flow are very similar to the constants of equation (l) 
which were evaluated for laminar flow. 

Curves of skin-friction ratio as a function of Mach number were 
obtained from equation (2) for the temperat\ire conditions of the data of 
figure 7 and are compared to these data in figure 10. The method of 
obtaining skin-friction ratio from equation (2) is explained in Appen- 
dix C. 

The agreement between the predicted values of skin-friction ratio 
from equation (2) and the experimentally determined values of skin- 
friction ratio is excellent over the entire range of Mach numbers and 
wall -temperature conditions for the experimental data of figure 7. It 
should be emphasized that only two experimental values were used to 
determine the constants, those at M^ = 3*82. The agreement at the other 
test conditions therefore represents a test of the method. 

Evaluation of skin friction for flight conditions .- The procedure 
described in Appendix C was used to obtain the curves of figure 11 from 
equation (2) where skin-friction ratio, Cp/Cp^, is plotted as a function of 
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Mach number, Mi, over a range of values of wall -temperature ratio, IV/Ti. 
The value of Ti = 392° R, corresponding to the standard isothermal 
altitude range of 35 >000 to 105,000 feet, was used for all calculations. 

The Sutherland viscosity law was used, with Si = 199° R, as given by 
the National Bureau of Standards (ref. I8) . A reference Reynolds nimi- 
ber, Ri = 10^, was used. 

There is a significant difference in skin-friction ratio between 
the zero -heat -transfer curves of figures 10 and 11 at Mach numbers 
greater than k. The zero -heat -transfer cirrve of figure 10 was calcu- 
lated for wind-tunnel-test conditions assuming a constant wind-tunnel 
reservoir temperature, therefore, Ti reduced rapidly at high Mach num- 
bers. The curves of figure 11 are more nearly representative of flight 
conditions at high altitudes. 

Comparison of the T* method with theory and experiment at subsonic 
speeds.- The curves of figure 11 predict veicy large changes in skin- 
friction ratio with changes in wall -temperature ratio at the low Mach 
numbers, particularly at subsonic speeds. Because of the lack of suit- 
able experimental data on flat plates with large rates of heat transfer, 
the predictions of the T' method were compared with experimental data 
on turbulent flow in smooth pipes at subsonic speeds. Since experimental 
skin-friction results at large rates of heat transfer in smooth pipes 
were very meager, heat-transfer data were also used in the form of 
Stanton number ratios, St/St/„ _ m ^ been demonstrated by 

Colburn (ref. I9) that Stanton number is proportional to skin friction 
or that Stanton number ratio is equivalent to skin-friction ratio. In 
figure 12, skin-friction ratio as determined by use of equation (2) at 
M = 0 is compared to experimental and theoretical results of skin- 
friction ratio and Stanton number ratio as a function of wall -temperature 
ratio for turbulent flow in smooth pipes . The Reynolds number based on 
pipe diameter, Rp, of the experimental data is 10^. Reynolds n\miber 
based on length of run of turbulent flow on a flat plate , Ri , was calcu- 
lated to be 2.6x10° for Rp = 10^, for a boundary-layer thickness equal 
to the radius of the pipe. This Reynolds number, Ri = 2.6x10® was used 
in equation (2) to obtain results for comparison with the data. 

Colburn (ref. I9) and McAdams (ref. 20) have shown that if viscosity 
is evaluated at a film temperature (equivalent in purpose to the present 
T»), Tf = 1/2 Tw + 1/2 Ti, the skin -friction equations for zero heat trans- 
fer in smooth pipes could be applied to flow with heat transfer. The 
Drew, Koo, and McAdams equation (ref. 19), 

Cp = 0.0014 + 0.125 Ri)-°-°^ 

and the Karman-Nikuradse equation (ref. 21), 

= 2 log3^o(RpvT4ci-)-0.8 
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vere ‘both evaluated at ttie film temperature^ as defined alove. The Drew^ 
Koo, and McAdams equation, and the Kannan-Nikuradse equation gave results 
at Rp = 10^ that differed very slightly from the results of equation (2). 
Skin-friction ratio, as determined from equation (2) agreed well with the 
experimental skin-friction results of reference 21 and was slightly higher 
over the entire range of wall -temperature ratios when compared to the 
heat -transfer results of references 21 and 22. It appears prohahle that 
equation (2) can he used to predict skin-friction ratio at subsonic Mach 
numbers . 

Comparison of the T* method with experimental results of refer- 
ence 23 »- For Mach numbers above 4*5 there is very little experimental 
data available against which to check the T* method other than the data 
already presented in the preceding sections of this paper. In fact, the 
only additional data in this Mach number range known to the authors are 
the data of reference 23. Unfortunately, these data are not for flat- 
plate conditions, being a set of measurements of the turbulent boundary 
layer on a nozzle wall at Mach numbers of 5-0^ 6.8, and 7- 7* The values 
of local skin-friction coefficients and the corresponding values of local 
incompressible skin-friction coefficients were based, in reference 23^ 
on the same (Reynolds number based on momentiim thickness) . It was 
therefore necessary for comparison with the present results to re-evaluate 
the local incompressible skin-friction coefficients so that the values of 
local skin-friction coefficients and the corresponding values of local 
incompressible skin-friction coefficients, Cf^, were based on the same 

Ri (Reynolds number based on length of run of turbulent flow on a flat 
plate) . The resulting expression for local incompressible skin-friction 
coefficient based on R^ and derived from the Karman-Schoenherr equa- 
tion is 


Cfn 


0,0293 


log 


2Rc 


10 


(cf/cf^) 


R. 


■a 


logio 


2Rc 


(cf/cfPj^ 



where Rg and (cf/cf.) 

^0 


are the values given in reference 23. 


The experimental results in the ratio of skin friction to incompres- 
sible skin friction, £>•^6 shown in figure I3 and are compared 

with skin-friction ratio as determined by use of equation (2) for the 
temperature conditions of the experiment. It can he seen that although 
the data show no definite trends of the effect of heat transfer on skin- 
friction ratio, the level of the data is in fair agreement with the 
results from the T’ method. The T' method predicts a large change 
in skin-friction ratio between zero -heat -transfer conditions and condi- 
tions where the wall -temperature ratio is unity; however, for the heat- 
transfer conditions of the experiment, the T* method predicts a rela- 
tively small change in skin-friction ratio. The measured changes in 
skin-friction ratio are, in general, of the same order of magnitude to 
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be expected from the T* method; however, in some cases, the changes 
are of opposite sign. 

Predicted effect of Reynolds nvimber on skin-friction ratio .- The 
predicted effect of Reynolds number on skin-friction ratio was investi- 
gated by evaluating equation (2) at standard isothermal altitude condi- 
tions, and using the procedure described previously. The effect of 
Reynolds nvimber on skin-friction ratio at Mach numbers of 2 and 8, at 
Tv = Ti and Tv = Tr, is shown in figure l4. It can be seen that at a 
Mach number of 2, over the wide Reynolds number range from 10 to 10 , 
that the predicted Reynolds number effect on skin-friction ratio is 
small. In fact, the magnitude of this predicted effect over the Reynolds 
number range from 3x10® to 30 x 10 ® is so small, less than 2 percent, that 
it would be difficult to confirm experimentally. This prediction is 
consistent with the results of Chapman and Kester (ref. l) where no 
Reynolds number effect on skin-friction ratio was detected at Mach num- 
bers up to 3.6. At a Mach nvimber of 8, however, the predicted effect 
of Reynolds number on skin-friction ratio over the Reynolds number 
range from 3x10® to 30x10® is the order of 5 to 10 percent. Over the 
Reynolds number range from 10® to 10®, the predicted effect is the order 
of 10 to 25 percent. This analysis was made to caution the reader about 
the validity of the assumption that skin-friction ratio is invariant with 
Reynolds number. Since the Reynolds numbers encountered in high-speed, 
high-altitude flight can be expected to be in the order of several hun- 
dred millions, the application of equation (2) with the assumption of no 
Reynolds number effect on skin-friction ratio may produce serious dis- 
crepancies . 


CONCLUSIONS 


Experimental measurements of skin friction of the turbulent boundary 
layer have been made on free -flying models at high supersonic speeds. 

The results of this investigation can be summiarized as follows: 

1. The effect of wall -temperature ratio on skin-friction ratio is 
large, an increase of the order of 35 percent at Mach numbers of 2.8l 
and 3 - 82 , when free-flight data of this experiment are compared with 
zero -heat -transfer skin friction data. 

2. Although skin-friction measurements at conditions of zero heat 
transfer are not available at Mach nimibers as high as "J, the results of 
this investigation indicate that the same trend of increasing skin- 
friction ratio with decreasing wall -temperature ratio will persist at 
very high Mach numbers . 

3. The T’ method of Rubesin and Johnson for laminar boundary 
layers has been used to evaluate skin friction of the turbulent boundary 
layer. By use of slightly modified equations with experimentally 
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determined constants, computed values of skin-friction ratio agree well 
with measured values over a wide range of Mach nijmhers and wall- 
temperature ratios. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Dec. 8 , 19 ^k 


22 


MCA TW 3391 


APPEiroiX A 

ESTIMATION OF SURFACE TEMPERATURE 


In order to determine the heat-transfer conditions of the present 
investigation, it was necessary to evaluate the wall -temperature ratio, 
Tw/T'i. The determination of T^/Ti was a transient-condition problem 
due to the short flight time of the models, the order of 0.01 second. 
Neglecting axial heat flow, the solution of this problem is given in 
paragraph Ul, equation (2) of reference 2k and is shown below in the 
nomenclature of this paper. 



where 

H average heat-transfer coefficient 
km thermal conductivity of the model material 
t time 

T absolute temperature 

initial temperature of the model 
Tr recovery temperature 
T-v, temperature at y = yi 

y normal distance from the midplane of the model wall 


yi half -wall thickness 

a thermal diffusivity of the model material 


(values tabulated in Appendix IV, 


3n positive roots of P tan 3 - 
ref. 24) 

For the present purposes, the series converges with one term. 
addition, the values of p^^ small so tan ps= 3 and 3i = J k^ 
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therefore, at y = Yi 


equation (Al) simplifies to 


Tw - Tj 
Tr - Ti 


2 


2 + 


yiH 


e 


oH 

yiV 


(A2) 


The values of Tv/Ti listed in this report are the mean values over 
the portion of the test models for which skin friction was measured. The 
maximum variation of over this portion of the models was 8 F vhich 

resulted in a maximum variation of T^/Ti of only 0.02. Axial heat flov, 
which was not considered, would have a negligible effect on T-^- in this 
region. 
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APPENDIX B 

DETERMINATION OF TURBULENT ORIGIN 


It was assumed that in the presence of small boundary -layer trips, 
the effective point of origin of turbulent flow could be determined by 
the following method: 



Assume the momentum thickness of the turbulent boundary layer, 0 iji, 
is equal to the momentum thickness of the laminar boundary layer, 9 ^, 
at the transition point, d. 


= ®Td 


Substitute 


then 




(CFx)od 





(Bl) 
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Nov* for the present investigation, the experimentally determined Cp 
is based on the area included between the indicated points b and c . 
The incompressible skin-friction coefficient based on the same area is 


CFi 




Ob 


(x) 


be 


so 


cf 

/be 


(Cpx) 


be 




(B2) 


bb 


If it is assumed that Cy/Cp^ is invariant with Reynolds number, then 
(cf/Cf^^ = (^Cf/Cf^^ ^ and equations (Bl) and (B2) can "be combined. 


(Cpx) 


Be 


Cfj^x 
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Cf^x 
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Cf^x 


oB 


Cf^x 


(B3) 


bd 


For each test condition, the position o can be found which will satisfy 
this equality. If Mi and Tw/Ti are known, the value of can Be 

oBtained from reference 25; it the Reynolds numBer per inch, Ri/x, is 
known, the values of can Be calculated from the Karman-Schoenherr 

equation. The value of Cp was determined experimentally. The left 

side of equation (B 3 ) , Tcp/Cp^^ , -vdiich includes the measured value of 

Cp, is not strongly dependent on the position of the effective turBulent 

origin; consequently, fCp/Cp. j as a function of the position o is 
nearly a horizontal line. The' right side of equation (B3)j, however, is 

very strongly dependent on the position o and intersects f Cp/Cp. j 

^ ^^Bc 

almost perpendicularly. The position of the effective turBulent origin 
is therefore sharply defined. 

The assumption that Cp/Cp^ is invariant with Reynolds numBer seems 
to Be justified Because existing data indicate that Cp/Cp^ is very 
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nearly independent of Reynolds number (ref. l) . According to the T* 
method of this report, it appears that there may be some dependence of 
Cp/CPi on Reynolds number; however, for the present purposes, the effect 
of Reynolds number on Cp/Cp^ would be small (the order of 2 percent) 
because of the limited range of Reynolds number variation. 


In the presence of large boundary-layer trips, the assumption of 
equal laminar and turbulent momentiim thicknesses at the transition point 
is unjustified; therefore, the following method was used to determine 
the position of the effective turbulent origin: 




Leading edge 




Effective turbulent origin 


/ 


Transition point 


Base of test model 



Base of tare model 


\ 




a o d 


At the transition pointy the momentum thickness of the turbulent 
boundary layer is equal to the momentum thickness of the laminar bound- 
ary layer plus an increment due to the presence of the trip. 


eTd=|(Cpx) 


od 


■ - ' V 


+ A0 


(Bl^) 


ad 


If A0 is related to the drag of the trip, it is proportional to only 
that part of the trip drag which results in removing momentum from the 
boundary layer. The term "trip drag” used hereafter refers to the per- 
centage of trip drag which results in removing momentiim from the bound- 
ary layer. The moment\mi loss in the boundary layer up to the point d 
is equal to the friction drag of the laminar boimdary layer plus the 
trip drag. 


Cp^ q 2itrx + q 2jtrx = J 2jtrpu(ui - u)dy 

o 


1 

2 








(B5) 


0 


NACA TW 3391 


27 


Combining equations (bU) and (B5) produces 



(b6) 


To relate this to some known quantity, such as the total drag, the trip 
drag is expressed as a percentage of total drag 


Cdt^^ = A(Cdx) 


ac 


(B7) 


Substitute equations (b 6) and (B7) into (b4) 




Cf/Cf^^Cf^x 


od^ 


Cf/Cf^ 


(Cf X ) + A(C])x) 

V ^ Ad 


ac 


bd 


Cf-x 


od 


(b8) 


Assume Cp/Cp^ is invariant with Reynolds number and combine equa- 
tions (B8) and (B2) . 


(Cfx) 


Idc 


CFl^^ ^ A(Cfx) 


ac 


Cf.x 


oc 


- ( Cf.x 


fcp X 


Ob 


V 


(B9) 


od 


To detemine the position o which will satisfy this equality^ it 
is necessary to make some assiimption regarding A. For the present 
investigation, it was assumed that all the momentum change due to the 
drag of the trip was confined to the houndary layer, and the trip drag 
was calculated hy subtracting the estimated components of drag from the 
measured total drag of the test models. Note that when the trip drag 
is negligible, A = 0 and equation (B9) reduces to equation (B3) • The 
same method was employed to solve equation (B9) as was used for 
equation (B 3 ) . 
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APPENDIX C 

PROCEDURE FOR DETERMINING SKIN-FRICTION RATTO 
BY THE T‘ METHOD 


In order to evaluate skin-friction ratio, Cp/Cp^, for any Mach 
number. Mi, at any wall -temperature ratio, T-w/Ti^ for a given Reynolds 
nijmber, Ri, by the T* method using equation (2) and the Karman- 
Schoenherr incompressible flow equation, the following procedure should 
be followed: 

mi , , 

1. Evaluate — from equation (2;, 

Ti 

1 ^ = 1 + 0.035 Mi^ + 0.45 - 

2. Evaluate R' from the following relationship: 




P’ Ti 

but considering constant pressure through the boundary layer, — = 

Pi ^ 

therefore , 



The ratio, (li’/pi), can be determined from the Sutherland equation. 


Pi ^ + sA 

Pi \Ti/ \T* + Si/ 

(Si = 216° R and Reynolds number, Ri = lo"^, were used to evaluate the 
curves of fig. lO) . 

3. Evaluate Cp' from the Karman-Schoenherr equation in the form 


0.242 

using R’ from equation (Cl). 


logio(Cp*R*) 


(C2) 
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4. Evaluate Cp from the relationship: 


Cf» _ Pi 
Cp P* 

Pi T' 

but since 

P* Ti 


Cp» _ T' 


(C3) 


using Cp' from equation (C2) and 
previously. 


Ti 


from equation (2) as determined 


5. Evaluate Cp^ from the Karman-Schoenherr equation in the form 


O.2U2 


logic (CF iRi) 


(Cl^) 


It should be pointed out that the evaluation of skin-friction ratio 
by use of equation (2) can be simplified if the Prandtl-Schlicting rela- 
tionship, 


Cp = 


0.46 

(logioR)^’® 


(C5) 


for incompressible flow is used instead of the Karman-Schoenherr equa- 
tion. 


When equation (C5) is combined with equation (2), the resulting 
equation for the evaluation of Cp/Cp^ is: 


Cp 

Cpi 


1 


logioRi 


logic 

(T'/Ti)(p»/ni) 


.6 


(c6) 


where Cp/Cpj^ can he evaluated directly. 

The authors chose to use the Karman-Schoenherr equation throughout 
this paper, although the Prandtl-Schlicting equation fits the existing 
low-speed skin-friction data equally well. The authors wished to he 
consistent with other authors who have used the Karman-Schoenherr 
equation. 
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TABLE I.- TEST CONDITIONS AND MODEL CONFIGUEATIONS 


(1) 

(2) 

( 3 ) 

( 4 ) 

(5) 

(6) 

( 7 ) 

(8) 

( 9 ) 

(10) 

Mo 

— 

X 

per in. 

Ml 

per in. 

Gun Lore 

i, 

in. 

Figure number 
for nose profile 
and boundary 
layer trip 

Tw 

Ti 

Tr 

Ti 

Symbol 

2 . 81 ^ 

1.61 

2.81 

1.58 

Rifled 

2.0 

2(a) 

1.03 

2 . 4 o 

0 

3.88 

2 . 24 

3.82 

2.13 

Rifled 

2.0 

2(a) 

1.05 

3.60 

□ 

7.09 

3.64 

^.63 

1.83 

Rifled 

2.5 

2(b) 

1.29 

6.64 

0 

7.31 

2.52 

6.90 

2.08 

Rifled 

2.0 

2(c) 

1.70 

9.47 

□ 

7 .i^l 

3-82 

7.00 

3.14 

Rifled 

2.0 

2(c) 

1.75 

9.72 

■ 

3.86 

2.23 

3.78 

2.13 

Rifled 

2.0 

2(d) 

1.05 

3.54 

D> 

3.80 

2.23 

3.67 

1.99 

SmootL 

2.0 

2 (e) 

1.05 

3 . 4 o 

<3 



TABLE II.- FINAL RESULTS 


Ml 

Tw 

Ti 

(0 

0 

R^lO-s 

RclO-® 

Cf 

Cp 

Symbol 

(see 

table I) 

2.81 

1.03 

1.58 

0.63 

3.00 

0.00284 

0.867 

0 

3.82 

1.05 

2.13 

.88 

4.07 

.00227 

.730 

□ 

5.63 

1.29 

1.83 

1.05 

4.71 

.00170 

.562 

0 

6.90 

1.70 

2.08 

1.03 

to 

2.97 

4.06 

to 

6.09 

.00125 

.404 

to 

.451 

D 

7.00 

1.75 

3.14 

1.35 

to 

5.21 

6.06 

to 

9.92 

.00115 

.395 

to 

.446 

1 

3.78 

1.05 

2.13 

1.75 

4.94 

.00204 

.694 

D> 

3.67 

1.05 

1.99 

.20 

3.78 

.00240 

.724 

<1 
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Figure 1 .- Typical models 
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(a) 10° half- angle wedge, 6 threads 0003 inch deep 



(b) Circular arc with 20° tangent angle, 7 threads 0.006 inch deep 




-0.12 


-0.10 • 


0.015 ? 


(e) 15° half-angle wedge, 6 threads 0.003 inch deep 


Figure 2.- 


Nose profiles and boundary- layer trips. 
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Figure 3»- Test model assembly, exploded view. 



36 


NACA TN 3391 



(a) 10° half -angle wedge, 6 threads O.OO 3 inch deep 



A-19634 

(h) 10° half -angle wedge, 6 threads 0.001 Inch deep 


Figure 1)-.- Shadowgraphs of test models at Mq = 3-90. 
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Mach number 
(a) Data of Reference 12 



0 12 3 4 

h Ri 

(b) Cross plot of 5(a) 

Rgure 5 The effect of Mach number on base 
pressure for two-dimensional wings in turbulent 
flow. 


Skin friction 

Incompressible skin friction 


UO 

00 



Mach number, M| 

Figure 6.- Effect of Mach number on skin -friction ratio with various wall- temperature 
ratios, with present data uncorrected for boundary- layer thickening at the trip. 
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o 
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cn 



Mach number, M| 


Figure 7.- Effect of Mach number on skin -friction ratio with various wall-temperature 
ratios, with present data corrected for boundary-layer thickening at the trip. 
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c 

o 

+- 
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c 

cn 



Tw/T, 



■P' 

o 


(a) M, = 2.81 

Rgure 8.- Comparison of experimental and theoretical results of the effect of wall- 
temperature ratio on skin-friction ratio at constant Mach numbers. 
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0 2 4 6 8 10 

Tw/T, 

(c) M, = 5.63 


Figure 8 



(d) M, = 7.00 


Concluded. 
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Skin friction 

Incompressible skin friction 
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Mach number, M| 


Figure 9 Comparison of skin-fricton ratio as determined 
by use of equation (1) with experimental values of skin- 
friction ratio. 


Skin friction 

Incompressible skin friction 



Mach number, M| 

Figure 10- Comparison of skin-friction ratio as determined by use of equation (2) with 

experimental values of skin -friction ratio. 
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Figure 11 - Effect of Mach number on skin-friction ratio at standard isothermal 
altitude, as determined by use of equation (2) at constant values of^ wall- 
tempera+ure ratio. 
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Figure 12.- Comparison of skin-friction ratio as determined by use of equation 
(2) with experimental and theoretical results at subsonic speeds. 
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SKin friction 

Incompressible skin friction 


-F 

On 



Mach number, M, 

Figure 13.- Comparison of skin-friction ratio as determined by use of equation (2) with 

experimental results of reference 23. 
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Skin friction 

Incompressible skin friction 
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(a) M, = 2.00 



I 2 3 5 10 20 30 50 lOO 


Reynolds number, R| , x 10‘® 
(b) M, = 8.00 



Figure 14.- Effect of Reynolds number on skin-friction ratio at standard 
isothermal altitude, as predicted by use of equation (2). 
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